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Abstract 

Charmed meson production in semi-inclusive deep inelastic scattering is investigated in the color 
dipole formalism. The transverse momentum distributions are calculated. We find good agreement 
with the HI data using a hard fragmentation function. 
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I. INTRODUCTION 



One of the most important heavy-flavor production processes is deep inelastic scattering 
(DIS). A quantitative understanding of heavy-flavor physics in ep and pp-collisions is a 
prerequisite for the discovery of new effects in pA and /L4-collisions. The color dipole model 
represents a good phenomenological approach to heavy-flavor production. In this model, the 
DIS cross section is factorized into a light-cone wavefunction, which describes the splitting of 
the virtual photon into a colorless quark-antiquark dipole with transverse separation r and 
a dipole cross section which depends on r and describes subsequent dipole scattering off the 
proton. The dipoles are eigenstates of the interaction and thus, multiple scattering effects 
on nuclei are easily described Historically, this was the initial intention to develop the 
dipole formalism. The dipole approach can be formulated in the target rest frame. For low-x 
the typical propagation length of a dipole fluctuation exceeds the interaction time by orders 
of magnitude. Because of time dilatation we can think of a frozen transverse separation 
during the interaction with the target. It is expected that the growth of the gluon density 
is slowed down at very low Xb by gluon-gluon recombination. Within the dipole approach, 
this saturation effect is described most naturally by introducing an x-dependent saturation 
scale. A fast-moving dipole produced from the incident virtual photon decouples from soft 
QCD interactions. This color transparency manifests itself in the dependence of the dipole 
cross section proportional to r 2 for small dipoles. The color dipole approach is inherently 
nonperturbative and includes higher twist effects which are important at low x-Bjorken and 
for low transverse momenta of the produced D-mesons. 

The paper is organized as follows. Section 2 deals with D-meson production in deep 
inelastic scattering. We develop the kinematic framework for the description of semi-inclusive 
deep inelastic scattering. Section 3 introduces the dipole approach to DIS and presents the 
virtual photon-proton cross section preparing the numerical evaluation of the deep inelastic 
cross sections. In section 4 we finalize the dipole calculation and compare with data from 
HERA's HI collaboration. Finally, in section 5 we summarize the results. 
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II. SEMI-INCLUSIVE D-MESON PRODUCTION IN DEEP INELASTIC SCAT- 
TERING 



In inclusive Deep Inelastic Scattering (DIS) we consider the process IN — > I'X, where 
a lepton / with momentum k scatters off the nucleon N with momentum P resulting in a 
momentum k' for the scattered lepton V . The target breaks up into an unobserved final 
state X. The cross section is most conveniently expressed in terms of the Lorentz invariants 

Q 2 = -q 2 = -(k-kr>0 v = ^ (1) 



or the ratios 



Q 2 P-Q , v 

XB= 2Mv V= p—k- (2) 



In the target rest frame, i.e. in the rest frame of the nucleon, v gives the energy transfer 
from the lepton to the nucleon and y the ratio of energy transfer to incident lepton energy. 
Q 2 is the virtual mass squared (virtuality) of the exchanged photon. We also use the center- 
of-mass energy squared in the / N system s and in the 7*p system 

S = (A; + P) 2 , W 2 = (q + P) 2 . (3) 

In Semi-inclusive Deep Inelastic Scattering (SIDIS) of the type IE -> VhX the 

hadron h is detected in coincidence with the scattered electron e' in the final state. The 
corresponding situation is pictured in Fig. [H The hadronic target H has mass M and 
four-momentum P, the produced hadron h has mass m and four-momentum ph. 

The SIDIS cross section can be expressed in the same Lorentz invariants xb and y as for 
inclusive DIS plus the final state hadron variables ph- Since we aim for pt factorization we 
split the final state hadron variables into the Lorentz invariant z momentum fraction of the 
hadron relative to the struck quark 

and the transverse hadron momenta ph ± relative to the 7* momentum direction. Both are 
invariant under boosts along the virtual photon axis. The differentials are related by 



FIG. 1: Kinematics of one-particle semi-inclusive scattering. The reaction plane in (e,e'h)- 
scattering is denned by the initial and final momenta of the lepton (k,k'). The produced hadron 
Ph and photon momenta q define a second plane. 



In analogy to inclusive DIS the SIDIS cross section can be obtained from the contraction of 
the hadronic and the leptonic tensor 

ZEhda _ 7ra 2 2My 

4 ^liv vv 



d 3 p h dx B dy Q 
2E " d0 (6) 



d 3 p h dndE' Q 

The conversion between these two forms is given by the Jacobian 

E' . 
dxndy = - —. — r=—dE d cos 6 

for azimuthal symmetry. L^ v is the spin averaged leptonic tensor and is the hadronic 
tensor 

x (2x)'6'(P + q-p,,-Px)- (7) 

The semi-inclusive hadronic tensor is related to the inclusive one via 

<n h (P,q) > W^(P,q) = J ^W, v (P,q, Ph ) (8) 

where < rih(P,q) > is the average number of particles produced of type h in a reaction 
defined by P and q. The above notation underlines the dependence on the momenta of 
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the deep inelastic reaction. Integrating over the hadron momentum ph one recovers the 
inclusive hadronic tensor times the multiplicity as shown on the right hand side of Eq. 8. 
Summation over all possible hadrons h would lead to the inclusive hadronic tensor. The 
tensorial structure of W^ u is determined by Lorentz and gauge invariance. 
Assuming time reversal and parity invariance we represent the electromagnetic hadronic 
tensor with four structure functions Wi(x, Q 2 , £,pf ) 

W» v {P, PhA ) = (^-^1 + ^ 

+ P ^ Tv + T " pL W 3 + ?k^W, (9) 
Mmh mr h 

where 

T M = — ^V- (10) 

The first two structure functions have the same origin as in the inclusive structure functions 
which are related to longitudinal and transverse photon scattering. The dependence on the 
hadronic angle 4>h allows two more structure functions. Next we introduce four dimensionless 
semi-inclusive structure functions H^x, Q 2 , z,p^ ± ). A factor of 2z arises in formal analogy 
to the known inclusive structure functions Fi 2 

2 

2zH 1 = M{W 1 + f^W 4 ) 



2zH 2 = HW 2 + §^W 4 ). (11) 



The integration over the azimuthal angle (ph removes the H 3 and structure functions 
from the semi-inclusive cross section 

Hi(x B ,Q 2 ,z) = \J dPh±_ d( i ) hH l {xBiQ 2 ,z ) p 2 hL ) 

= 7T J dp^HiixB^^Z^U (12) 



which leads to 

da 87ra 2 ME 



x B y 2 H 1 (x B ,Q 2 ,z,p 2 h± ) + (1 - y)H 2 (x B ,Q 2 ,z,p 2 h± ) 



(13) 



dx B dzdy Q A 

In the high-energy approximation E = s/(2M) the cross section has a form similar to the 
inclusive cross section. Factorization of deep inelastic scattering and hadronization allows to 
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represent the inclusive structure function H 2 as a product of the known quark distribution 
functions and the quark fragmentation function D^/gJz) jj| 

H 2 (x B ,Q 2 ,z) = H 2 (x B ,z) = ^ e 2 f [qf/ H (x B )D h / qf (z) + qf /H{x B )D h / qf (z)}. (14) 

/ 

This factorization has been proven to leading order in Q p|. Here, in this paper, we want to 
calculate the structure functions Hi in the dipole formalism for 7*p. Using equation (jSJ), we 
connect the Hi to the standard structure functions for inclusive lepton-hadron scattering, 

< n h (x B , Q 2 ) > F x {x B , Q 2 ) = J dzH x {x B ,z,Q 2 ) 

< n h {x B ,Q 2 ) > F 2 (x B ,Q 2 ) = J dzH 2 (x B ,z,Q 2 ). (15) 

Following the inclusive case, we need two structure functions specified for longitudinal and 
transverse photon polarization. Polarization interference in semi-inclusive DIS occurs only 
for W3 and W4 which vanish after integration over the azimuthal angle. 

W L = £ (\ = 0)-W-e(\ = 0) = 2z(J^H 2 -^j (16) 

1 2z 
W T = -[e(X = -l)-W-e(X = -l)+e(X = +l)-W-e(X = +l)} = —H 1 

where e is the virtual photon polarization vector. For the virtual photon momentum q = 
(q°, 0, 0, q z ) we define the polarization vectors: 

e(X = 0) = i (q\ 0, 0, q°) e(X = ±1) = -L(0, 1, ±i, 0). (17) 

The combined longitudinal and transverse electron-proton cross section is 

da na 2 M ( f y 2 



dxB dz iy dpi± - 2QV I V 2 + 1 " ») Wt + f 1 " ^ )■ < 18) 
We recall the relation between the total 7*p cross section and the inclusive hadronic tensor 
W^(P,q) 

<4^x(A) = ^^e;(X)W^(P, q)e u (X). (19) 

Since the amount of D-meson production is measured by the branching ratio 

1 f , da™^ DX {X) 



<n D >= tot / dy ' (20) 

<*V.x(A) J dy 

we obtain 

2E D da^ DX 8vr 2 aM , /xXTrruiy/ „ N 

T3 = < A )W» V {P,q,p D )e v A , 21 

d 6 p D ys * 
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or with 



ysz 



da L > T 



-L,T 



ii 2 Ma em d 2 pr> dz 
We insert this expression in (1T81) to find 



W ■£ 



L.T 



(22) 



a,. 



do[ep -> DX) 
dx dy dz d 2 p^ nxy 

+ (1- 



2 



-y+1 



da T (Yp -> M) 



2/) 



da L (Yp^DX) 



(23) 



This equation relates the SIDIS cross section to the semi-inclusive cross section of a 
virtual photon on the proton. The latter can easily be calculated in the color dipole model, 
since the photon is the source of the color dipole, and the target proton mainly delivers the 
low x-gluon which is necessary for the virtual heavy quark-antiquark state in the photon to 
materialize. 



III. D-MESON PRODUCTION IN THE COLOR-DIPOLE MODEL 



The 7*p charm quark production cross section is calculated in the dipole model (cf. ref. 



2211) . We start with the photon- gluon-fusion processes shown in fig. EJ The respective 



Feynman graphs are represented by the amplitudes Mf and 





b) 




V(k T )Y° 




V(k T )Y° 



FIG. 2: Boson-Gluon- Fusion graphs. FIG (2a) shows iMi, FIG (2b) iM 2 . 



iMS = E^AK^eeqT* ^ (24) 
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Here e(X) is the polarization vector of the virtual photon 7* with polarization A, u aq (p q ) 
the spinor of the heavy quark with helicity a q and momentum p q , and v a Apq) the spinor 
of the antiquark with helicity a q and momentum p q . The charge of the heavy quark in 
units of e is eg. The polarization vectors are given in eqs. (1T71) . The target color field 
represented by the vertex function V{k t ) depends on the transverse momentum k t injected 
by the gluon from the target. The vertex V is restricted to its time-component, i.e. the 
formula holds only in the target rest frame. It is assumed that the target gluon mediates 
transverse momentum only, longitudinal momentum and energy transfers are negligible. In 
the high-energy approximation the fermion propagators can be treated onshell [5J]. This 
approximation is equivalent to neglecting the instantaneous terms in light-cone quantization 
6j. We can express the energy denominators in terms of LC variables. The coordinate 
system is chosen such that the z-axis points in (f-direction. With the light-cone momentum 
of the photon (q + ,0±,q~) = (Q 2 / {2mx b) , 0, — ttixb) and a as the light-cone momentum 
fraction of the quark we find (p q ,p q _L,P q ) = (o!q + ,p q j_, {p 2 q ± + m 2 )/(aq + )). This allows to 
simplify the propagators in eq. (24). Equivalent expressions hold for the antiquark with a 
replaced by 1 — a 

1 —a 



(q - p q ) 2 - m 2 p 2 ± + e 2 
1 a — 



(q - p q ) 2 - m 2 (p q± - k t ) 2 + e 2 

e 2 = a(l-a)Q 2 + m 2 . (25) 

The parameter e 2 is related to the inverse extension squared of the quark-antiquark state. 

The Yp cross section is proportional to the absolute square of the summed transition 
amplitudes Mf . In the phase space integration we use the integration variable k t given by 
the equation k t = p q + p q — q instead of p q \ 

dot d 2[ v^~ d 2 ki 

d"o= £ ; 8(27T)5(g o)2 a(1 _ ^ (AMA'XMf + M-)(Mr + M 2 -.) (26) 

The polarization-interference terms vanish when integrated over the azimuthal angle. 
It is advantageous to transform all expressions into impact parameter space since at high 
energies the quarks propagate along fixed impact parameter trajectories. We express the 
vector amplitudes M^{k u p q ) by their Fourier transforms M? (6, r) in eq. (27) 
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M?(kt,p q ) = J d%dh-M?(b,r)e i >* tUi ^. (27) 
The explicit form of the vector amplitudes in coordinate space is: 

Mf(6,f) = -i2 Pq >ee Q (l - a) J d ^0^e-^-^V(k t ) 

(p q ± - h) 2 + e 2 

With the vertex function in impact parameter space 



and the light-cone wave function: 



*»M = e QV ^VKV^a) I ^ e- ^WT^r^ (29) 
one finds that the momentum space integrals factorize and yield the convenient expres- 



sions: 



M?(b,r) = - t 2p° q ^^V(b + r)^(r) 

(30) 

M^lr) = i2p°^^V(b)W(r). 
The light-cone wave function describes the splitting 7* — > qq and appears identically 
in both initial state amplitudes. The vertex function has the respective arguments b and 
b + r which are the impact parameters of the antiquark in M 2 and M 1 . Using p° q ~ aq°, 
p° q ~ (1 — a)q°, we find for the 7*p production cross section of a heavy quark with momentum 

Pq 



da 



L.T 



J _J_ J d 2 bd 2 r 1 d 2 r 2 dae i ^ n - r ^ (31) 



x 

CqCq 



9 



Due to the integration over the antiquark momentum which is equivalent to the k t integra- 
tion, the impact parameters in (Mf + M^) and (M™ + M™) coincide. Since the momentum 
of the quark p q remains unintegrated, the exponentials of the Fourier transforms survive. 
The light-cone wave functions simplify after summation over spins and polarizations 

(32) 

^(r 1 )^ T (f 2 ) = iE^E^^WjfWr^). 

Combining the gluon vertex operators V appropriately we can construct dipole cross 
sections 

d 2 b[V(b) - V{b + n)][V\b) - V\b + f 2 )} = 

~ J d 2 b [| V(b)-V(b + n) | 2 + | V(b)-V(b + f 2 ) | 2 - | V0+r f 2 )-V(jb + f f 1 ) | 2 ] .(33) 

The corresponding absolute squares can be grouped into dipole cross sections employing 
the definition: 



j d 2 b\V{b)-V{b + r)\ 2 . (34) 



One finds the following final form for the differential 7*p cross section for charm quark 
production, i.e. q = c: 



da ' (p/*p -> (cX)) 1 



d? P i (2tt; 



J d 2 n d 2 r 2 da e^*-**) V L < T (n)y* L > T \ 



r 2 



x ^{cr c5 (fi) + a C c{r 2 ) - a c£ (fi - r 2 )}. (35) 

A similar expression has been derived in [3] to describe the transverse momenta of Drell- 
Yan pairs. In the differential cross section light-cone wave functions enter with Bessel func- 
tions Kq and K\ depending on e = \/m 2 + a(l — a)Q 2 and r% or r 2 : 

2N a e 2 

^ L (a,n)^ L (a,r 2 ) = ° ™ Q AQ 2 a 2 (l - a) 2 K (e ri )K (er 2 ) 

( 27r ) 

V T (a, n )** T (a,r 2 ) = Q -^^{ m 2 Q K Q {er{)K»{er 2 ) 

+ e 2 [a 2 + (1 - ayp^K^K^)}. (36) 

T\T 2 
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If one integrated out the transverse quark momentum in equation (1351) . one would find 
a delta function in (7*1 — r 2 ) and recover the total cross section of charm production. The 
differential cross section does not have such an intuitive interpretation as the total DIS 
cross section. As shown in the derivation a physical quark-antiquark fluctuation is not a 
prerequisite for the appearance of the dipole cross section in the inclusive cross section. 
Instead, interference terms in the square of the transition amplitude produce dipole cross 
sections when a quark in the amplitude and a quark in the complex conjugate amplitude 
enter with different impact parameters. This becomes very obvious in the Drell-Yan process, 
where the dipole formula for the transverse momentum distribution of lepton pairs looks 
very similar to (1351 ) although there is not necessarily any dipole fluctuation involved in the 
Feynman graphs [8|. 

We gave the inclusive DIS cross section in terms of the cross section of a quark-antiquark 
pair scattering off a proton. For the calculations in this paper we use a fit to the dipole cross 
section and its energy dependence. For an explicit derivation of the dipole cross section 
from QCD we refer to the loop-loop correlation model, see Ref. Q|. In the following we use 
the GBW-saturation model by Golec-Biernat and Wusthoff [l(3] with a dipole-nucleon cross 



section of the form: 

cr w( a; > r ) = °o 

where 



6XP I I I (37) 



< 38) 

This cross section is not a well-defined quantity for Q 2 = 0. Therefore, x has been modified 

/ Am 2 n \ 

x -> x y 1 + -qt ) ( 39 ) 

in order to describe the transition to the photoproduction region. The distance between 
the quark and antiquark is r. The inverse momentum 1/Qq = 0.2 fm is a perturbative 
distance. At x = xq = 0.41 * 10~ 4 twice this distance equals the distance where the cross 
section changes from the perturbative r 2 behavior into a constant behavior. The parameters 
A, do and xq are given in Table (JIJ) [nil]. The fit is performed within a Bjorken interval 
4 * 10~ 4 <x B < 0.01 and virtuality range of O.lGeV 2 < Q 2 < 400GeV 2 . 

The saturation model contains four fit parameters. This is a small number compared to 
the number of parameters in parton distribution functions. The authors have performed a 
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cr (mb) 


A 


x 


l/Qo(/m) 


Oqq with charm 


29.12 


0.277 


0.41 • 10~ 4 


0.2 



TABLE I: Fit parameters in the GBW saturation model [10]. The charm quark mass is fixed to 
m = 1.5GeV. 



fit to HERA data and reached a satisfactory x 2 in Ref. [10j containing charm in the sum 
over flavors. 

The model describes two different saturation effects. The GBW cross section saturates 
for quark-antiquark pairs with large r separation. Arbitrarily large gg-separations are not 
physical due to hadronization, but this does not matter, since they occur in the qq light- 
cone wave function with very low probability. The other saturation effect is implemented 
via the x-dependence of Rq ~ l/Q s > where Q s acts as a saturation scale. The 7*p cross 
section rises earlier to a for decreasing x. One expects this increase to be slowed down or 
stopped at even lower x when the gluon-gluon-recombination cross section becomes sizeable 



111 ]. The color-glass model predicts a new QCD domain which is nonperturbative but not 
dominated by confinement effects. Note, both saturation effects are not to be mixed up with 

n 

the saturation of the gluon density in impact parameter space discussed in Ref. [9|. 
For our calculations we can expand the cross section to first order in r 2 



<7qq{X, r) ^ 



<7Q 



ARl{x) 



cr r 



(40) 



where R is given in (j38i) . In the following we refer to this use of the color dipole cross 
section as r 2 -approximation. For D-meson production with a charm mass m = 1.5 GeV, we 
do not expect that the size of the semi-inclusive cross section is very much influenced by 
large dipole separations. The dipoles contributing have a size of about 0.2 fm. 



IV. NUMERICAL CALCULATION 



With the analytical formulae and the above dipole cross section we can numerically 
calculate the double differential charmed quark production cross section da T ' L /d 2 p^r. We 
reduce the fourfold Fourier integral over r*i, r 2 which occurs in the 7*p sub cross section ( !35l) 
to a one-dimensional integral over the dipole separation r. We find integrals with modified 
Bessel functions of the second kind K 01 and Bessel functions of the first kind J 0) i 
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with 



h = J drr J (p^r) K (er) a q g(r) 

J drr 2 J (pjr) K^er) a q g(r) 

h = J drr.h{p^r) /^(er) <r M -(r) . (42) 

The small r 2 -approximation ( 1401 of the GBW dipole cross section allows to perform the 
r-integrations J1-/3 analytically 



h . 
J2L 



do 4(e 2 -p, 



<r 16e(e 2 -2pf) 

-approx. 4jR 2 ^±2 + £ 2)4 



3 1 r 2 -approx. ^j£2 fo±2 _j_ £ 2^3 ' V ' 

The integrals require an extension parameter e greater than zero to be finite. For heavy-quark 
production with e > m this is no problem. We have checked how much the approximate 
recalculation and the full GBW dipole cross section differ for the sum of transversely and 
longitudinally polarized photons. In Fig. [3] we show the j*P -* cl cross section in r 2 - 
approximation over the same cross section calculated with the full GBW parametrization. 
Deviations appear mainly at low transverse momenta of less than about 1 GeV 2 . We find 
good agreement for higher Q 2 , when the average dipole size is small. Above p^r 2 > 2 GeV 2 
we find a very small deviation of maximally 2%, which is sufficient for our calculations. 

Next, we study the relative contributions of longitudinally and transversely polarized 
virtual photons in Fig. 01 The transverse cross section is clearly dominating. The relative 
contribution of the longitudinal cross section grows when the virtuality of the photon changes 
from Q 2 = 2 GeV 2 to 100 GeV 2 . This behavior is similar to vector meson production. 
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0.96 1 1 1 1 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 18 20 

p t 2 (c) [GeV 2 ] 

FIG. 3: The ratio of the cross section da{^*p — > cX)/dpjr 2 in r 2 -approximation over the same cross 
section calculated with the full GBW cross section a{r) is shown as a function of the transverse 
momentum of the charm quark p^r 2 . 



To calculate D-meson cross sections we must let the charm quark fragment. The fragmen- 
tation function Dj*(z*) gives the probability that the original charm quark with a momentum 
P fragments into a D-meson with momentum fraction z*P. All our momentum fractions in 
section 3 refer to the photon momentum q. Since we start with a charm quark of momentum 
P = aq and end up with a D-meson with momentum zq the momentum fractions multiply: 



z = z a. 



(44) 



To calculate D-meson production we convolute the charm quark production cross section 
da L ' T ( r y*p — > cX) with the nonperturbative fragmentation function D(z*): da L,T ( , y*p — > 
DX) 



da 



L.T 



dz d 2 p D 



i^p^DX) 



dp^: da da 



L.T 



a d 2 Pc da 

x ir D (-) 5U--P. 

\aJ V a 



( 7 *p - cX) 



We use the Peterson fragmentation function 



12j which has the form 



(45) 
(46) 



D h Q (z*) 



n(h) 



z*[l - 4 



£ Q 12 
l-z*J 



(47) 
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FIG. 4: The charm quark production cross section 7*p — > cX induced by transverse and longitudinal 
photons for Q 2 = 2 GeV 2 and W 2 = 200 GeV 2 . The transverse photon cross section is dominating 
over the whole transverse momentum region. 



This fragmentation function has been successful in describing experimental data of heavy 
flavored hadrons and can be easily understood in the following way. The transition amplitude 
is proportional to the inverse energy difference AE between the initial heavy quark Q and 
the final meson h plus light quark state q 

AE = Eq — E h — E q 



m 2 +P 2 - x m 2 + z 2 P 2 - Jml - (1 - z) 2 P 2 



m 2 a / I 



where = uiq can be used for simplicity. The fragmentation function is given by the square 
of the transition amplitude normalized correctly. The Peterson parameter cq = m 2 q jm 2 Q 
parametrizes the hardness of the fragmentation and is fitted to data. The larger the quark 
mass, the smaller is eg and the harder is the fragmentation, i.e. Dq(z*) is peaked more 
towards z* = 1. The z-integrals of the Peterson fragmentation functions Dq are normalized 
to the branching ratios fq which summed over all hadron species add up to unity 

£# = £/ D h Q (z*)dz* = 1. (49) 



h h 

In general we are using the branching ratios given in Ref. 



13l | which refer to the specific 
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kinematical domain of the Hl-experiment. The world averages evaluated by the Particle 
Data Group lj| are also listed in Table (JTTJ) for comparison. The large isospin violating 
branching ratio for the Z)°-meson compared with the D + -meson branching ratio comes from 
the accidental nearness of the D*-resonances to the sum of the masses of the D-meson and a 
pion. Namely both the D* + and D*° can decay into D° with large branching fractions, but 
the D* + can only decay to the D + with a small branching. The Peterson parameter for the 
charm fragmentation function is e c = 0.05 in the leading-logarithmic approximation (LLA) 



15j. Although the branching ratios into D-mesons and _D*-mesons are approximately equal, 
it is much easier to measure D*- mesons than D-mesons at HERA, since _D*-mesons have 
a good signature due to their D* — > Dn — > Kixix decay chain. The fits are performed with 
data from 'clean' e + e~ reactions. We use these functions in DIS assuming independence 
from the underlying production process of the heavy quark. 

D° D + D*° D*+ 



PDG 0.565 ± 0.032 0.246 ± 0.020 0.213 ± 0.024 0.224 ± 0.028 
HI 0.658 ± 0.054 0.202 ± 0.020 0.263 ± 0.019 

TABLE II: Branching ratios c — > D for several -D-mesons. Reference values from Particle Data 
Group PDG Hi and HI Q. The HI data is fitted without theoretical constraints due to the 
experimentally covered kinematical region. 



The fragmentation process of bottom quarks requires a harder fragmentation function 
since the b-quark mass is 4-4.4 GeV, i.e. about 3.5 times larger than the charm quark 
mass. The reference value e& = 0.006 for the bottom quark is therefore much smaller in the 
leading logarithmic approximation. The attachment of a light quark degrades the bottom 
quark momentum less than the charm quark momentum. The LLA-Peterson fragmentation 
functions are plotted for both flavors in Fig. O 

When a high-momentum charm quark in the photon fragments into a charmed meson 
without additional transverse momentum, the scaling variable z* gives the transverse and 
longitudinal fractions p^r /Pd = Pc IPd = V- 2 * a t the- same time. One finds: 

- 1 dz* 



a mean ratio of transverse momenta < p^ IPd >= 1-68 and < p^/Pb >= 1-29 for the bottom 
sector, respectively. This makes the differential cross sections of the heavy quark harder than 
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* 

z 



FIG. 5: The Peterson fragmentation function with e c = 0.05 for charm quarks and = 0.006 for 
bottom quarks in the leading log approximation. The bottom quark fragmentation is harder since 
the bottom quark mass is about 3.5 times higher than the charm quark mass. Both curves are 
normalized to unity. 



the meson cross sections by the corresponding ratios. In order to see how the fragmentation 
acts we plot the longitudinal Yp cross section da L /d 2 p 2 for the two different values of €q 
as a function of the transverse momentum in Fig. El One sees how harder fragmentation 
functions shift the curve to higher transverse momenta. 



V. COMPARISON WITH EXPERIMENTAL DATA 

We evaluate numerically the ep cross section ( 1231) for the HI experiments 0, Q at 
HERA taking into account the experimental cuts. In order to specify these cuts we use two 
different coordinate systems, the K-system and the K*-system: The z-axis of the K-system 
points in the incoming electron beam direction. The z-axis of the K*-system points in the 
photon direction. The experimental results of the measurements by the HI collaboration 
are given in the K-system, i.e. with respect to the electron beam. The photon virtuality 
varies 2 < Q 2 < 100 GeV 2 , the inelasticity 0.05 < y < 0.7 and pseudorapidity \r)(D)\ < 1.5. 
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FIG. 6: The longitudinal heavy meson production cross section j*p — > h(Q)X induced by longi- 
tudinal photons as a function of the transverse momentum of the meson. The two e parameters 
correspond to charm (e c = 0.05) and bottom quark (e& = 0.006) fragmentation. For comparison we 
also show the transverse quark production cross section without fragmentation. 



Further the minimum transverse momentum pjj min is 2.5 GeV for the center-of-mass energy 
=319 GeV 13] and p^ min = 1.5 GeV for y/s = 300 GeV H- The theoretical dipole 



cross sections we use include well this Q 2 -range. Our calculation, however, extrapolates 
somewhat the Bjorken-x range used in the GBW dipole cross section. The GBW range 
includes the x^-interval [10~ 4 , 10~ 2 ], the experimental data spread over the larger interval 
[10 _5 ,0.02]. Therefore the theoretical results for the smallest xb values should be regarded 
with caution. We convert the differential cross section for electron-proton scattering (l23l to 
the experimentally given variables. For x,y,Q 2 we use the relation s = m 2 + Mp + Q 2 / (xy), 
i.e. Q 2 = xys in the high-energy approximation. The rapidity is approximated by the 
pseudo-rapidity, y rap , ~ 77. In the target rest frame (TRF) where p% = ttit cosh rj and 
v = ys/(2M), the z-fraction of the D-meson is calculated as follows 

P-Pd trf P° D 2M 

z = — = — = rriT cosh 77 . (50) 

P ■ q v ys 

We evaluate the differential -D + -cross sections with respect to the photon virtuality 
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FIG. 7: Theoretical cross section da(ep — ► D + X)/dQ 2 with branching ratio fj? + = 0.202 plotted 
against experimental data from reference [13]. Note, the experimental lower cut for the transverse 



momentum is pj=, , 



2.5 GeV. 



da/dQ 2 and the rapidity of the D + -meson da/drj. For both cross sections we integrate 
over the p^ range inside the experimental cuts. In Fig. we show da(ep — > D + X)/dQ 2 



for e c = 0.05 and m c = 1.5 GeV and compare with the HI data [I3j. The Q 2 -decrease is 
determined by the integrals (c.f. eq. (42)) over the Bessel functions K 0> i(er) depending on 
the inverse extension parameter e 2 = a(l — a)Q 2 + m 2 . With increasing photon virtuality 
Q 2 , the Bessel functions damp the size of the cross section. 

In Fig. [8] we plot the da(ep — > D + X)/drj cross section against the HI data [131 ] . For the 
total cross section, higher quark masses would result in a lower rj distribution. This holds 
for the whole rapidity range, even more so at mid rapidity. The whole curve becomes higher 
for a harder fragmentation. This is expected since the differential cross section falls off for 



high rapidities of either sign. Better data are announced [17j to come out for D*-mesons. In 
order to obtain theoretical predictions for other D-mesons one just has to multiply the D + 
cross section with the appropriate ratio of branching ratios (c.f. Table II). For example by 
multiplying the D + cross section with 2.6 one obtains the combined D* + , D*~ cross sections. 

Figure [9] shows the theoretical da(ep — ► D*X)/dpj^ cross section against the HI data fl^ ]. 
Our theoretical cross section reproduces the data adequately. It tends, however, to under- 
shoot the cross section at large px- This would mean that the gluons in the proton have too 
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FIG. 8: Theoretical cross section da(ep — > D + X)/drj with branching ratio / ( 



against experimental data from Ref. 
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small transverse momenta or the unintegrated gluon density Jg{x, k±) underestimates large 
transverse momenta. Indeed, the GBW parametrization for the dipole cross section leads to 
a Gaussian fc^-dependence for the unintegrated gluon density, while a power dependence is 
expected in QCD. 

We have also investigated the dependence on the hardness of the fragmentation process. 
Our calculation is not very sensitive to changes of the Peterson parameter e c . A harder 
fragmentation lifts the differential cross section for transverse momenta above ~ 1.5 GeV 
and reduces it for lower ones. 



VI. SUMMARY AND OUTLOOK 



In this paper the color dipole model has been applied to the production of heavy flavored 
mesons in semi-inclusive deep inelastic scattering. In particular, transverse momentum dis- 
tributions are a more valuable source of information about the production and fragmentation 
processes than total charm production. In the target rest frame, the virtual photon splits 
into a heavy quark-antiquark pair which subsequently scatters off the target. A color dipole 
is the lowest Fock state of the virtual photon. Instead of considering higher Fock states con- 
taining gluons, which become more important for lower x, these contributions are absorbed 



20 




0.01 ' 1 1 1 1 1 

2 4 6 8 10 



PtP*) [GeV] 

FIG. 9: Theoretical cross section da(ep — > D* + X)/dpj^ with branching ratio ff?* + = 0.263 plotted 
against experimental data from reference [16l |. 

into the x-dependence of the dipole fit. The dipole formulation is in a mixed representa- 
tion, where the longitudinal momentum is treated in momentum space and the transverse 
position in coordinate space. For high energies, the quark-antiquark pairs with fixed sepa- 
ration are interaction eigenstates. This formulation is very convenient for treating multiple 
scattering effects. With the help of the framework developed in this paper one can also 
calculate nuclear effects on heavy quark production by replacing the dipole-proton cross 
section by the dipole-nucleus cross section in Glauber approximation. For electroproduction 
on the proton we find a differential cross section in Q 2 which agrees well with HI data. 
The large phase space coverage helps to collect sufficient statistics in the experiment. On 
the other side it presents a real challenge for the theoretical parametrization of the dipole 
cross section. In this respect the parametrization of Golec-Biernat and Wusthoff's has a 
clear advantage compared with other approaches due to its wide kinematical applicability. 
It would be interesting to see the effects of the improved GBW dipole cross section j3| on 
the calculations. Inclusion of these corrections into the dipole description of direct photon 
production [lj|] reduces the cross section at high-p^. Nevertheless, our results suggest that 
the dipole cross section is not the main source of deviations. 

The fragmentation function has a stronger influence on the shape of the transverse mo- 
mentum distribution. We get good fits with a rather hard fragmentation function. Since 
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there is no previous work which combines hadronization with the dipole approach, it is 
worthwhile to investigate how to treat heavy quark production and fragmentation consis- 
tently. Our fitted dipole cross section contains higher twist effects, i.e. it includes many soft 
gluons in the partonic cross section. Therefore we propose to use a fragmentation function 
at a low scale without evolution. We do not consider evolution which may become important 
at p^ > m c . 

Intrinsic transverse momenta of the target gluons are encoded into the dipole cross section. 
A primordial transverse momentum of the projectile gluon is often introduced in parton 
model, even in NLO calculations to harden the spectrum. Usually the primordial momentum 
substantially exceeds the typical hadronic scales. This is not necessary in the color dipole 
model, which generates the intrinsic gluon transverse momenta automatically (e.g. see in 



Similar calculation can be applied to hadron induced production of charm 



211]. While the 
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total cross section of charm production in pp collisions is well explained 
seen how well pt distribution can be explained in the same framework. Once it is possible 
to reduce the errors arising from hadronization, one may obtain more detailed information 
about the accuracy of the underlying dipole calculation. 

Within the same dipole formalism one can calculate charm production in diffractive DIS, 
I + p — > V + cc + p. This cross section turns out to be a higher twist, ~ l/ m c- However, 
semi-inclusive diffractive production, I + p — > I' + ccX + p, is a leading twist, ~ V m c 
due to additional gluon radiation. It was demonstrated recently that diffractive hadro- 
production of heavy flavors is also a leading twist, and dipole model calculations well explain 
available data. 
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APPENDIX A: EXPERIMENTAL CUTS 



The experimental cuts are given with respect to the electron beam line. We relate the 
transverse momenta relative to the electron beam line and the virtual photon. As shown in 




FIG. 10: Relating transverse momenta relative to the electron beam line, pj^ , and to the virtual 
photon direction pjj for a small rotation of the z axis. 



Fig. [10] we can write 



-»-]_' ->± , z -< 
Pd = Pd + -It 
a 




and therefore 



(Al) 



Pd = y q't) + (Ph) 2 ~ 2 Z ~ qiph cos0 , (A2) 
where the transverse virtual photon momentum q' t is given by kinematics 



, Q 2 M 

% = — \iQ^( 1 -y)- 1 ( A3 ) 



~ Q^l^j. (A4) 

The momenta are shown in Fig. [TU1 The quantities relative to the electron direction are 
indicated with a prime. The transverse photon momentum is accompanied by a factor 
z/a = z* . This factor evolves because the virtual photon splits into a heavy quark pair and 
not directly into D mesons. 
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